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Abstract 

We present an inflationary scenario based on a phenomenologically viable model with 
direct gauge mediation of low-scale supersymmetry breaking. Inflation can occur in the 
supersymmetry-breaking hidden sector. Although the reheating temperature from the in- 
flaton decay is so high that the gravitino problem seems to be severe, late time entropy 
production from the decay of the pseudomoduli field associated with the supersymmetry 
breaking can dilute gravitinos sufficiently. We show that gravitinos are also produced from 
the pseudomoduli decay and there is a model parameter space where gravitinos can be the 
dark matter in the present universe. 
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1 Introduction 



Cosmological inflation [TJ [2] in the early universe is now considered as a part of the "standard" 
cosmology. It can solve many cosmological problems such as the horizon problem and the 
flatness problem. It also accounts for the origin of the primordial fluctuations [3] . In most 
models of inflation, a scalar field with the flat potential (called the inflaton) which drives 
inflation is assumed to be outside the standard model of particle physics since there is almost 
no inflaton candidate in this framework^] Then, a natural question is how inflation dynamics 
is embedded in a particle physics model beyond the standard model. To construct an inflation 
model, it is reasonable to respect supersymmetry (SUSY), which is one of the most promising 
candidates for physics beyond the standard model, because the radiative corrections for scalar 
fields are suppressed in supersymmetric theories and an inflaton potential can be naturally 
flattened. 

If SUSY is realized in nature, it must be broken at some energy scale. We usually leave 
its dynamics to a SUSY breaking hidden sector and the SUSY breaking effect is transmit- 
ted to the visible sector by some interactions. Gauge mediation is an attractive mechanism, 
which uses the standard model gauge interactions [3, E\ (For recent general arguments, see 
Ref. [9]). It naturally suppresses the unwanted flavor-changing processes due to flavor blind- 
ness of the gauge interactions. While many gauge mediation models are already known, the 
model proposed by Kitano, Ooguri and Ookouchi (KOO) [10J is a distinguished one (see also 
Ref. [IT]). This model is based on the recent development of metastable SUSY breaking in 
supersymmetric QCD (SQCD) started from the work of Intriligator, Seiberg and Shih (ISS) 
[T2] and is classified as so-called direct gauge mediation models [13] where the flavor symme- 
tries of a hidden SUSY breaking sector are weakly gauged and identified with the standard 
model gauge symmetries. The KOO model can generate sizable gaugino masses while many 
direct gauge mediation models suffer from anomalously small gaugino masses [13] (see also 
Ref. |15J). However, since the vacuum considered in the KOO model (the ISS vacuum) is 
not the global minimum of the potential, it is necessary to impose several conditions on the 
model parameters in order to guarantee the vacuum stability. In particular, two hierarchical 
mass scales are required in the flavors of this massive SQCD model. We must also consider 
the cosmic history that selects the ISS vacuum in this modelj^l 

*It is usually considered that Higgs field cannot be responsible for inflation since its potential is too steep 
to account for the primordial fluctuations [4]. Recently, however, several inflationary models that regard the 
standard model Higgs boson as the inflaton were proposed [5]. In these models, the Higgs self-coupling is 
effectively suppressed by increasing the effective Planck mass or the Higgs kinetic term and hence primordial 
fluctuations consistent with the observational data [5] can be generated. 

2 One approach is to consider the finite temperature effect on the SUSY breaking sector [16]. For instance, it 
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Recently, two of the present authors realized hybrid inflation [T7] embedded in a gauge- 
mediated SUSY-breaking model motivated by the KOO model [18]. This scenario uses two 
hierarchical mass scales in the magnetic description of the KOO model such that the higher 
scale corresponds to the inflationary scale and the lower scale corresponds to the SUSY 
breaking scale. Since the inflaton rolls down to the ISS vacuum after inflation, it naturally 
explains why the SUSY breaking vacuum is selected in the cosmic history. Another attractive 
feature of this model is that if we can find the model parameter space that accommodates 
the overall cosmic history (including inflation and the present dark matter abundance) and 
gauge mediation, we might be able to prove (or disprove) the scenario by the near future 
collider experiments and the cosmological observations such as Large Hadron Collider (LHC) 
and PLANCK pTO]. 

In this paper, we present an inflationary scenario following the line of Ref. [18] and analyze 
the reheating stage of an inflation model with gauge mediation. Since the couplings of the 
inflaton in the SUSY breaking hidden sector with the visible sector fields are given by the 
standard model gauge interactions, the reheating process after inflation is predictable in this 
framework. We here take a simpler Wess-Zumino model than that in Ref. [18] as the SUSY 
breaking hidden sector, but the essential feature is the same. As is the case in Ref. [20] (see 
also Ref. [21]), the evolution of the universe after inflation proceeds as follows. After the 
end of (hybrid) inflation, the damped oscillation of the inflaton and waterfall fields dominate 
the energy density of the universe and the Hubble parameter of the universe decreases. At 
the stage of reheating from the inflaton decay through the messenger loops, only the fields 
in the visible sector are thermalized, since those in the SUSY breaking sector are heavy. 
As the Hubble parameter decreases, the pseudomoduli associated with the SUSY breaking 
that was stabilized at the origin of the field space during inflation starts to oscillate around 
the metastable vacuum f| Because the lifetime of the pseudomoduli field is rather long, the 
pseudomoduli oscillation dominates the universe. Then, the pseudomoduli decays before the 
Big Bang Nucleosynthesis (BBN). Gravitinos produced in the thermal bath at the reheating 
stage are diluted enough by the entropy production from its decay, whose present abundance 
otherwise tends to be so large that it overdoses the universe. Although gravitinos are also 
produced in this decay process, the branching fraction is suppressed since the pseudomoduli 
can decay into the visible sector particles through the messenger loops. Thermally and non- 
thermally produced gravitinos can both contribute the dark matter abundance in the present 



is known that the metastable vacuum is preferred than the SUSY preserving vacuum at the high temperature 
in the ISS model. However, this approach has a difficulty that gravitinos are overproduced when the SUSY 
breaking sector is thermalized. 

3 It depends on the model parameters that the moduli oscillation takes place earlier or later than reheating. 
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Table 1: 



The charge assignments of the fields in our model. 



universe. 



The rest of the paper is organized as follows. In section 2, we present our inflation model 
and analyze the vacuum structure and the mass spectrum around the metastable vacuum. 
We also evaluate the soft mass spectrum of the visible sector fields from gauge mediation. In 
section 3, we present the inflationary scenario based on Ref. [IB] and estimate the reheating 
temperature after inflation. In section 4, the oscillation of the pseudomoduli field is analyzed. 
Next, we investigate the decay temperature of the pseudomoduli field. We will also show 
the gravitino abundance due to the pseudomoduli decay. Then, we find the model parameter 
space consistent with the BBN constraints and the present dark matter abundance. In section 
6, we conclude our discussions. 



In this section, we present a metastable SUSY breaking model, which drives inflation and 
gauge mediation, and analyze its vacuum structure. Then, we analyze the mass spectrum 
around the SUSY breaking vacuum. We also show the soft mass spectrum of the visible 
sector fields obtained by gauge mediation. 

2.1 The model 

The model is a Wess-Zumino model with an SU (N) global symmetry^ The matter content 
is summarized in Table [TJ \ and x are singlets under the SU(N) group and the fields p, p 
belong to the (anti-) fundamental representation under the SU(N) group. Z and Z are also 
a vector-like pair of the (anti-)fundamental representation under the SU(N) group. Y and 
$ are singlets under all the global symmetries. The tree-level Kahler potential of all fields is 

4 Later, we gauge this global symmetry to be the standard model gauge symmetry. 



2 SUSY breaking 
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assumed to be canonical. The superpotential is given by 

W = m 2 Y + p 2 ^ - h YX Y X - h^p^p - h z ( X Zp + pZ X ) - m z ZZ, (2.1) 

where the mass scale m is assumed to be larger than the scale p, and hy, h$, h% are coupling 
constants. The model has two 17(1) global symmetries. U(1)r symmetry is explicitly broken 
by the last term in the superpotential, whose breaking size is determined by the mass scale 
m%. First, we treat them as free parameters and will determine their sizes later as they pass 
all the conditions. For simplicity, we assume that all the couplings and the mass parameters 
are realB 

Let us analyze the vacuum structure of our model. There exists a metastable SUSY 
breaking vacuu nfl where the expectation values of the fields are denoted aqj 

— 777 

Y = p = p = Z = Z = 0, x = x =—=, (2.2) 

y iIy 

and the singlet <3> is the pseudomoduli of the SUSY breaking vacuum that is massless at the 
tree-level, but obtains a nonzero mass from radiative corrections. Since the vacuum energy is 
given by Vo = p A at the tree-level, the gravitino mass is estimated as 

m3 " = vlb (2 ' 3) 

where Mpi ~ 2.43 x 10 18 GeV is the reduced Planck mass. Conversely, when we express the 
SUSY breaking order parameter p in terms of the gravitino mass, it is given by 

,^ 7 . 9xl0 »GeVx( I ^) 1/2 , (2.4) 

where we take 15 GeV for the gravitino mass as the reference value. 

The pseudomoduli $ is stabilized by the 1-loop effects of the massive modes p, p, Z and 
Z. The supersymmetric mass terms for these massive modes are given by the following mass 
matrix: 

(p,z)m(£) = {p,z)(~'£* 




Integrating out these modes, we can derive the 1-loop effective Kahler potential for the pseu- 
domoduli field $ as [23] 



i-eff 



1 1 32tt 2 



MM* ( log 



(2.6) 



5 The sizes of the mass parameters can be explained by introducing additional gauge dynamics 
6 In the absence of the last term in (|2.1[) , SUSY is always broken in the model. 



7 In order to obtain the equal expectation values of the \, \ fields, we have used the additional condition 
from the D-term potential when we gauge the U{\)\ symmetry. 
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where A is a cut-off scale. Substituting the expression of the mass matrix (12.51) into this 
formula, we can obtain the following effective Kahler potential at the leading order of the 
scale m z : 

N 



i-eff 



32vr 



h*m z ($ + $ f ) +hl\<5>\ 



Note that the U(1)r breaking terms in the above expression are proportional to the mass 
parameter mz- The 1-loop effective scalar potential for the pseudomoduli field $ can be 
derived from the effective Kahler potential as 

2 



( d 2 K eS > 


-1 


dW 


V^$9$t ) 




<9$ 



(2.8) 



Then, we can extract the vacuum expectation value and the mass squared of the pseudomoduli 
field $ from the above effective scalar potential, which are given by 

. . 1 777-2 , 

' °'~27T' arg$0 = ' 

2 r N h Y h% // _ 2 ( 2 ' 9 ) 

o47r / h% m A 

As we will see, the field configurations finally settle down to this vacuum after inflation. Later, 
we consider the gauge mediation effects on the visible sector fields at this vacuum. 

There is a SUSY preserving vacuum far away from the origin of the pseudomoduli field 
space, in addition to the metastable vacuum. The vacuum expectation values of the fields are 
given by 

- m 2 _ /i 2 - h\ mV 

XX=j~, PP=J-, ZZ = h h m 2 » 

tiy n,0 nyriq, m z ,^ ^ 

hyh^, mz ' hyhq, mz 

The SUSY breaking vacuum has a nonzero transition probability to the SUSY preserving 
vacuum, which can destroy the stability of the SUSY breaking vacuum and hence SUSY 
breaking mechanism may not work. However, its decay rate T vac is evaluated by using the 
triangle approximation as [23] 

T vac oc e~ s , (2.11) 



h% 



hyh% \m z J V 
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Fermions 


Bosons 




Weyl mult. mass SU(N) 


Real mult. mass SU(N) 


$ 


1 1 


2 m cw 1 


Y, X, X 


1 0(y/h^m) 1 
1 0{y/k^m) 1 

1 9vik 1 


2 0{yfh^m) 1 
2 O(v^ym) 1 
2 ^ 1 


Z, Z, p, p 


27V 0(fem) □ + □ 
27V 0(^=m) □ + □ 


47V O(^m) □ + □ 
47V O(^m) □ + □ 



Table 2: The mass spectrum and the representations under the SU(N) symmetry. In the 
table, gy denotes the U(l)i gauge coupling. The scalar component of the pseudomoduli field 
$ has a mass (12. 9p by the 1-loop effects. 

Then, if the mass hierarchies m 3> \i and m ^> mz are realized, we have S 3> 1. Thus, 
the decay rate is sufficiently suppressed and hence the stability of the metastable vacuum is 
guaranteed. 

2.2 Mass spectrum 

We next analyze the mass spectrum around the SUSY breaking vacuum. The masses of the 
fields in the SUSY breaking sector except for the pseudomoduli are generated at the tree-level. 
The results of the masses as well as the representations of the fields under the symmetries are 
summarized in Table |5J The mass spectrum of the Y, Xi X sector fields are supersymmetric 
at the tree-level. All of the scalar fields in this sector except the Nambu-Goldstone mode 
have 0{y/hYm) masses. The corresponding fermionic modes have the same masses at the 
tree-level. The Nambu-Goldstone mode is associated with the spontaneous breaking of the 
U(l)i symmetry due to the nonzero expectation values of the Xi X fields. In order to give this 
mode a nonzero mass, we can weakly gauge the U (l)i symmetry. Then, the Nambu-Goldstone 
mode is absorbed into the massive gauge boson of the symmetry group. When we gauge the 
spontaneous breaking global symmetry, the fermionic partner of the Nambu-Goldstone boson 
becomes a part of the massive vector multiplet. 
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The Z, Z, p, p fields have the SUSY breaking mass spectrum and hence play the role 
of messengers in gauge mediation. The fermion mass matrix is given in ( 12. 5 p and they all 
have 0(^=m) supersymmetric masses which correspond to the messenger scale in our gauge 
mediation model. The SUSY breaking masses of the scalar components are small compared 
to the supersymmetric ones. 

As shown in the previous subsection, the scalar component of the singlet field $ obtains 
a nonzero mass due to the 1-loop effects. The fermion component corresponds to the would- 
be goldstino mode associated with the SUSY breaking, which is absorbed into the massive 
gravitino when the theory is promoted to supergravity. 

2.3 Gauge mediation 

We next consider the soft mass spectrum of the visible sector fields. In order to transmit 
the SUSY breaking of the metastable vacuum to the visible sector by direct gauge mediation, 
we embed the standard model gauge group into the SU(N) symmetry of the model. The 
formulae of the leading order gaugino and scalar masses are given by [25] 

2 Q 

mA -ii^ F *^ logdetM ' 

(2.13) 

,2 _ V- 



where gi (i = 1,2,3) is the gauge coupling constant of the U(l) x SU(2) x SU(3) group of 
the standard model gauge symmetry and the factor C l 2 is a quadratic Casimir. M is the 
mass matrix of the messenger fields p, p, Z, and Z which is given in (I2.5p . M s denotes the 
eigenvalue of this matrix. Then, at the lowest order of the scale mz, the gaugino and scalar 
masses in our model are explicitly given by 



g 2 h Y h^ u 2 m z 



167r 2 h\ mm 1 
/_Z ^ 2 U6W hi m 2 ' 



As we can see in the above expressions, the leading order gaugino masses do not vanish. Since 
the metastable vacuum is the higher-energy state in the potential, this result is consistent 
with the general theorem presented in [H]. Now, we parametrize the ratio between the gluino 
mass nig and the right-handed slepton mass m g for later use as follows: 

r g = nig/nig. (2-15) 
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In our scenario, this ratio is constrained by the requirements to accommodate the proper 
inflation, BBN and dark matter abundance in the present universe. 

We can rewrite the coupling h$ in terms of the observable parameters by using the ex- 
pressions of the moduli mass, the soft scalar mass, and the gaugino-to-scalar mass ratio r g as 
follows: 



We can see from this expression that the coupling h$ is an O(10 2 ) quantity when the 
observable mass parameters have the reference values. 

3 Inflationary scenario 

In this section, we study the hybrid inflation model embedded in the above theory. We see 
that the Y scalar field acts as the inflaton and the \i X fields act as the waterfall fields. We 
analyze its dynamics and density fluctuation generated during inflation. We also study the 
decay of the inflationary fields after inflation. 

3.1 Inflationary dynamics and cosmological perturbation 

First, we see the realization of hybrid inflation which can occur on the pseudo-flat direction 
(different from the pseudomoduli $ of the metastable vacuum) in our model. In order to do 
so, we concentrate on the following field configuration of the system: 

p = p = Z = Z = $ = 0, (3.1) 

which respects the global symmetry of the model. The relevant part of the superpotential 
(12. ip is, now, 

W ~ m 2 Y - hyxYx- (3.2) 

This type of superpotential is exactly of the form for F-term hybrid inflation [17] . That is, 
the field Y is the inflaton and the fields \i X are the waterfall fields. The tree-level scalar 
potential in the global SUSY limit is, then, given by 



tree 



hYXx\ +h\\Y\ 2 {\ X \ 2 + \x\ 2 ). (3.3) 



The masses of the scalar fields \i X are m s = ^yI^I 2 ^ hym 2 . For the large field value of Y, 
the fields Xi X are stabilized at the origin of their field space and hybrid inflation with the 
Hubble parameter, 

3^ < 3 - 4 > 
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-> X 
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Figure 1: The shape of the scalar potential V ra d is shown. It is convex upward and hence the 
slow-roll parameter rj will be negative. 



can occur. On the other hand, for the small field value of Y, \Y\ < Y c = m/V^Y, one of the 
scalar modes is destabilized. At that time, \ and x ac t as the waterfall fields and roll off to 
the point near the metastable vacuum discussed in the previous section with nonzero vacuum 
expectation values. 

Let us study the radiative correction, which drives the inflaton motion. Here, we redefine 

y 

Y = ~^J^- %6 an d consider the motion of the inflaton y. The potential for the phase 9 remains 
flat and hence it does not have any influence on the inflaton dynamics. The mass splittings 
between the bosonic and fermionic partners of x an d X due to the SUSY breaking result in 
the 1-loop effective potential for y as 



Kadfe/) 



+ 



327T 2 



+ h Y m 2 I log ( 1 



2m 2 
h Y y' 



hYm 2 I log ( 1 



2m 2 



(3.5) 



where M* is a cut-off scale. Defining a dimensionless parameter as x = y/y c 
(y c = \/2Y c ), the total potential is rewritten as 

V(<f>) ~ ^trcc + V iad 



h Y y/V2- 



m 



+ (1 



32tt 2 



2 log 



/ hyfn 



2 x 2 



X 



2\2 



log 1 + — 



V M 2 



1 



x 



2\2 



log 1 9 



(3.6) 



The shape of V ra ^ is shown in Figured! As long as h\/?>2Tx 2 C 1, we can approximate the 
potential energy with the tree-level value V ~ m 4 and the dynamics of y is governed by the 
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1-loop potential Kad- in this approximation, the slow-roll parameters are given by 



e = 



V 



v 



2567r 4 m 2 ' 



[x 2 + 1) log 1 



X 1 



[X 



l)log 1 



X* 



-i 2 



V" hi Ml 



t Y Jla Pl 



V 16ff 2 m 2 



x- 



(3x 2 + 1) log 1 + — + (3x 2 - 1) log 1 - — 



X/ 



(3.7) 



where the prime denotes the derivative with respect to y. Note that e is much smaller than 
\rj\ due to the suppression by a factor of hy/lQff 2 . The slow-roll conditions, e, |?7| <C 1, are 
satisfied until x ~ 1 -v^> y ~ y c , provided that hy 2 /2y/2ff <C m/Mp\. In this case, inflation 
ends when ?/ reaches to y c and the waterfall fields become tachyonic. Note that supergravity 
correction does not change the result significantly. Hereafter, we consider such a case. 

Before proceeding the discussions, we comment on the stability of this inflationary trajec- 
tory. In order for the above discussion to be valid, all the fields other than the inflaton must 
be heavy enough. Let us consider the scalar potential in supergravity, 



DiWK %3 DjW* 5-IWI 

3 Mp Y 



(3.8) 



dW/dfa + (dK/d^W/M^ and K** is the inverse 



where K is the Kahler potential, D$V 
matrix of d 2 K / dtfiidcf)*. During inflation, the term in parenthesis of Eq. (I3.8P is evaluated as 
~ 3H 2 Mp V As a result, the fields (collectively denoted ip) with the canonical Kahler potential 
receive the Hubble induced masses, 



(3.9) 



and this is true for almost all the fields in the theory. Thus, during inflation, such fields are 
well stabilized at the origin of their field space. In particular, the pseudomoduli field $ does 
not lie directly on the metastable vacuum, but lies at the origin during and just after inflation. 
This fact is not important for the inflationary dynamics but affects the cosmic history through 
the late time oscillation as we will see in the next section. 

We now evaluate the cosmological perturbation in this model. The number of e-folds M 
is calculated as 



1 f y V J 



16n 2 m 2 



x'dx' 



on m 
h\M 2 x 



1), 



(3.10) 



where we have approximated with log[l ± 1/x 2 ] ~ ±x 2 for larger field values^] From this 



This approximation is correct with an accuracy of 0(1) in the region in which we are interested. 
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expression, we obtain the field value evaluated M e-folds before the end of inflation, 

T 




—mxj^f ~ 
fly 




h 



m 



Y 



Z7T 



for M < 
for N > 



87r 2 m 2 

h^M 2 ' 
8ir 2 m 2 

h Y M 2 { 



(3.11) 



We find that a trans- Planckian initial value can be avoided if \FlJhym <C Mpi for the upper 
case and y/Jfhy < 2ir for the lower case. We do not need to worry that the inflaton may 
roll off to the SUSY preserving vacuum when we start from such a large initial value of the 
inflaton field, because the pseudomoduli field $ is trapped at the origin of the field space. 
The primordial density perturbation generated during inflation is, then, evaluated as 



V. 



1/2 



n 



(- 



H 



2e V2vrM P1 



4V67T 



m 



3 

2_ IN 
h Y V 3 

,1/2 



h% 2 Ml 



pi 

m 



Mi 



pi 



for A/" < 



for M > 



8ii 2 m 2 

8ir 2 m 2 
h\M 2 { 



(3.12) 



For the COBE/WMAP normalization, V]i ~ 4.9 x 1(T 5 at = 0.002Mpc _1 (6], we require 



m 



5.9 x 10 15 GeVx < 



h 



3 x 10- 3 

AcOBE 



51 



1/3 



-1/4 



for h Y < 3 x 10~ 3 , 
for /iy > 3 x 1(T 3 , 



(3.13) 



where Acobe is the number of e-folds after the COBE scale leaves the horizon. Note that 
when there is the moduli dominated era, Acobe is expressed as 



A^cobe ~ 51 + -lo^ 



5.4 x 10 14 GeV 



- log 

3 6 



10 10 GeV 



log 

3 S 



T, 



dom 



0.4 GeV") + 3 



log 



Ta 



4MeV 



(3.14) 



where Td om and are the cosmic temperature at the time of the moduli domination and that 
at the moduli decay, respectively. We will see the validity of the reference values in the next 
section. The spectral tilt can be evaluated as 



27r 2 ^n 2 



n s = 1 — 6e + 2rj 

Acobe 

and the scalar-to-tensor ratio r is given by 



1 

0.98, 



for 
for 



hy < 3 x 10~ 3 , 
h Y >3x 10~ 3 , 



(3.15) 



16e 
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10/3 
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16tH 

h\ 



4 /6 ^ P i 

m 



k 27T 2 A/cOBE 



for h Y < 3 x 10' 



for h Y > 3 x 10' 



(3.16) 
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which is much smaller than 0.1. Hereafter, we consider the case with hy < 3 x 10 3 and hence 
they have not been ruled out by current observation [6] 



3.2 Reheating after inflation 

We here investigate the reheating temperature due to the decay of the inflaton and the 
waterfall fields. Before proceeding with the discussion, we express the coupling hz in terms 
of the observable mass parameters for later use. By using (I2.4p . the moduli mass (12.91) . the 
expression (12 . 1 6[) of the coupling hq> and the estimation (13 . 1 3j) . it is given by 

-1/3 



1.8 x 10' 



h 



Y 



. (3.17) 



2 / rn 3/2 \ / mg \ 
v 3.5/ V 15 GeV/ V300GeV/ Vl.5TeV/ \3 x 10"^ 

Note that the coupling hz has an 0(1O~ 3 ) value when the observable parameters are set for 
the reference values. Then, the messenger scale can be estimated to be (9(10 13 ) GeV. We 
can also rewrite the mass parameter mz in terms of the observable parameters by using the 
gaugino mass formula fl2.14j) as 



m z ~ 8.2 x 10 12 GeV x 



rg\ 3 f m 3/2 \ 



771$ 



(3.18) 



fN\3.5J V 15 GeV/ V 300 GeV/ Vl.5TeV, 

which gives the mass scale of the explicit U(1)r symmetry breaking. 

After inflation, the oscillation of the waterfall fields X = \ + X an d the inflaton dominates 
the universe. Then, the decay of these fields reheats the universe. These fields have 0(^/hYm) 
masses as shown in Table |2] and (almost) maximally mix with each other to form the mass 
eigenstates. They dominantly decay into MSSM gaugino pairs through the messenger loop if 
the decay channels to the messengers are closed, 2hz > hy- The decay diagram is given in 
Figure [2j We can estimate the decay width as 

IN 4 /1d 9X 2 

r R 



amy 

2 (4vr)5 { hi J 



hY"m. 



(3.19) 



The reheating temperature is, then, estimated as 



90 



0.45 x 



1/4 



X v^rMf! 

1/2 



f h 



Y 



(4vr) 2 V 8tt 
5.2 x 10 10 GeV x N 7/2 



^1 (mMpi) 1 / 2 



hi 



3.5 



15 GeV ) 



300 GeV ) 



m~ g 



-V 



1.5 TeVJ \3xl0- 3 



h 



Y 



17/3 



(3.20) 



The tensor perturbation is hard to be detected even by the future observations. 
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Figure 2: The decay of the waterfall field X and the inflaton Y into an MSSM gaugino pair. 
■0 denotes the fermion component of the messenger multiplet. 

where gf- is the number of the relativistic degrees of freedom at the time of reheating and we 
take it as gf 1 ~ 220. The decay into the pseudomoduli field is also possible but is suppressed 
by the additional small Yukawa couplings. Since the messenger scale is estimated to be 
O(10 13 ) GeV, the SUSY breaking sector cannot be thermalized in our scenario. 

As well as the visible sector fields, gravitinos are also produced in the thermal bath at the 
reheating stage. From the above reheating temperature, the abundance of gravitinos is given 

by m 

(th) 

^^^^OeVx^H^)-'^), (, 21 , 

where s is the entropy density. Note that gravitinos are also produced by other processes 
such as the inflaton decay, however these are sub-dominant processes if the coupling hy is not 
so small so that we can ignore these contributions to the gravitino abundance. On the other 
hand, the present dark matter abundance is given by [6] 

n DM h 2 ~0.11, (3.22) 

where h = i?o/(100 km sec _1 Mpc _1 ) and Hq is the present Hubble parameter. Then, we have 
the constraint on the gravitino abundance, 

^<^1^4.1xl0- 10 GeV. (3.23) 

s s 

Since we would like to have the 0(1 — 10) GeV gravitino mass and the 0(1) TeV gluino mass 
in our scenario, we can see from Eqs. ( 13.211) that gravitinos are overproduced at the time of 
reheating after inflation. Then, some late time entropy production is required for this scenario 
to be successful. Actually, this model has a possible late time entropy production mechanism, 
that is, the decay of the pseudomoduli field. As we will see, it can dilute gravitinos sufficiently. 
Thus, we may obtain the correct abundance ( I3.22p of the gravitino dark matter. We study 
this issue in further details in the next section. 
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4 Gravitino dark matter 



In this section, we consider the cosmological evolution after reheating and investigate the 
gravitino abundance, which is similar to the case of Ref. [20]. First, we discuss the oscillation 
of the pseudomoduli field around the SUSY breaking metastable vacuum and the domination 
in a cosmological stage. Next, we will investigate the decay temperature of the pseudomoduli 
field. Then, we calculate the gravitino abundance due to the pseudomoduli decay. From the 
result of the decay temperature, we also estimate the amount of the entropy production by 
the moduli decay and evaluate the present abundance of thermally produced gravitinos. We 
find the model parameter space consistent with the constraints obtained from the analysis of 
these quantities. 

4.1 Moduli oscillation 

We here analyze the entropy production due to the moduli oscillation. As shown in the 
previous section, the moduli field is stabilized at the origin during inflation. As the Hubble 
parameter decreases after inflation, the potential minimum of the pseudomoduli moves away 
from the origin. Then, it starts to oscillate when the Hubble parameter becomes smaller than 
the moduli mass (12. 9p . The temperature at that time is given by 



/ go \V4 



(4.1) 

1.2 x 10 10 GeVx ( m * ^ 



300 GeVJ 

Here, g° sc is the number of the relativistic degrees of freedom at the onset of the pseudomoduli 
oscillation and we take it as g° sc ~ 220. Since the lifetime of the moduli field is rather long as 
we will see, the pseudomoduli oscillation dominates the energy density of the universe. The 
temperature when the domination of the moduli oscillation occurs is given by 



- dom 





f 


( |$ol 


s* dom , 








f 


( |*o| 


s dom , 







^osc for T R > T osc , 

(4-2) 

T R for T R < T osc , 



where g dom is the number of the relativistic degrees of freedom at the pseudomoduli dom- 
ination. When the pseudomoduli field decays, a large amount of entropy is produced and 
gravitinos are diluted. The dilution factor A" 1 due to the pseudomoduli decay is expressed 
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as 



1$ 



0| 



for T R > T osc , 
for T R < T osc , 



; t r Vv^Mp!, 

where Td is the decay temperature of the pseudomoduli, which we will estimate later. 

We here comment on the stability of the pseudomoduli oscillation. By using the observable 
mass parameters, the expectation value of the moduli field on the metastable vacuum is 
rewritten as 

The messengers become tachyonic when the moduli field takes the value, 



h Y h$ m z " V15GeV7 \1.5TeV, 

In order for the field configuration to settle down finally to the metastable vacuum, we require 
the following condition on the gaugino-to-scalar mass ratio: 

r g < 4.5. (4.6) 

Compared to the case in Ref. [20], the amplitude of the oscillation is rather small. Then, the 
abundance of gravitinos produced at the time of reheating is comparable to that of gravitinos 
due to the moduli decay. 

4.2 Moduli decay and DM abundance 

Next, we analyze the decay temperature of the pseudomoduli field and the abundance of the 
gravitino dark matter. The pseudomoduli can decay into both the visible-sector particles and 
gravitinos through the messenger loops which are encoded in the effective interactions of the 
pseudomoduli field with the visible-sector fields. We consider the case that the pseudomoduli 
field can decay into two Higgs bosons, i.e., m$ > 2rrih. We will see the case that the moduli 
mass is small and the field cannot decay into two Higgs bosons in Appendix. 

The interaction Lagrangian of the pseudomoduli field with the scalar components of the 
visible-sector superfields can be extracted from the moduli dependence of the soft scalar mass 
squared, m^($), which is given by 



9 2 i \ \h Y h% /i 4 , 3 h\h% fi 4 m z , + 



^ z V 16tt 

i 



h\ m 2 4 h\ m 4 



(4.7) 
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where the first term is the leading order of the scale and leads to the soft masses shown in 
section 2. The second term is the next-to- leading order one. There exist higher order terms 
as well, but we do not write them explicitly. The effective interactions of the moduli field can 
be extracted from this expression as 

dm 2 A<&) 

£ /=T* //t + hc - 

The pseudomoduli field dominantly decays into two Higgs bosons through the messenger loop. 
Other scalar fields are too heavy for the pseudomoduli field to decay into them. The decays 
into the standard model gauge bosons are sub-dominant because the effective interactions 
of the pseudomoduli field with the gauge bosons (shown in Appendix) are suppressed by a 
1-loop factor more than the interaction given above. The decay width is estimated from this 
effective interaction (14. 8[) as 

where in and Xg £1X6 defined as 

(iM + (iM • i - 6 * i(r6 < (4 - io) 

We have used the expressions of the observable mass parameters such as the gluino mass m g . 
The moduli decay produces large entropy in the universe and the temperature after the decay 
can be estimated from the decay width (14 .9p as 

T d ~ v /T u M Pl 

g \-2/ m g y ( m 3/2 \ -1 / ( 4 - 12 ) 
.3.5/ U.5TeV/ V15GeW VSOOGeVy 

The temperature is required to be above 2 MeV so that the standard BBN properly occursj^ 
This requirement constrains on the parameters in our model. 

We next calculate the number density of gravitinos from the pseudomoduli decay. The 
decay into the longitudinal mode of gravitino (would-be goldstino) is the dominant process. 



10 Baryogenesis must take place before the BBN. The required baryon-to-photon ratio is ns/n 7 ~ 6 x 10 -10 
[6]. In the case that the moduli oscillation dominates the energy density of the universe, larger baryon 
asymmetry must be generated in order to take into account the entropy production from the moduli decay. 
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The interaction Lagrangian with the longitudinal mode can be read off from the effective 
Kahler potential (12. 7p . 

£ ^-0^ir S 2 * , A/* /2 + h.c, (4.13) 



where ^3/2 denotes the gravitino. Then, the partial decay width of the pseudomoduli field 
can be calculated as 



1 / m 



2 \ 2 
<t> 



where we have used the expressions of the observable mass parameters. Compared to the total 
decay width ( 14.9ft . the decay width into gravitinos (14. 14[) is small by the quartic dependence 
of the coupling h®. The number density of gravitinos is expressed as 

n 3/2 3 T d 

— £- = - B3/2 x 2, 4.15 

s 4 m$ 

where B^^— ^3/2/^ h) is the branching fraction into two gravitinos and s is the entropy 
density of the universe. The branching fraction is given by 

,2 /„ \ 2 
18 (4tt) 4 x| 

Then, the density parameter of non-thermally produced gravitinos is estimated as 

tf* 0.033 x f^Vf-^V^f_^_r S . (4.17) 



5 3/2 - 7^ 77^ ^ ( — ) • (4.16) 



m 9 



3/2 « ^.u^xy ^GeVj 

This quantity is required to be below the present dark matter abundance, fi DM h 2 ~ 0.11. 
This gives another constraint on the parameters in our model. 

Since we have calculated the decay temperature of the moduli field, we can estimate the 
dilution factor A -1 due to the moduli domination before the decay. We here require that 
the reheating temperature is smaller than the temperature at the onset of the pseudomoduli 
oscillation, Tr < T osc . The reason is as follows. If the reheating temperature after inflation is 
too high, the abundance of thermally produced gravitinos is too large to explain the present 
dark matter abundance even if we include the effect of the dilution by the pseudomoduli 
decayl"] Then, from the condition Tr < T osc , the constraint on the coupling hy can be 
expressed in terms of the observable parameters as 

h < 22xl n-3 x 1 ( r § V 8/17 ( ^2 \W ( X 21/34, m v- 18/17 

11 If the reheating temperature is too low, we cannot neglect the direct gravitino production [27]. We do 
not consider such a case. 
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5 10 15 20 25 10 20 30 40 50 60 



m 3/2 /GeV m 3/2 /GeV 

Figure 3: The allowed region (shaded region) of the gravitino mass (horizontal axis) and the 
gluino mass (longitudinal axis). In the left and right figures, the moduli mass is taken to be 
300 GeV and 500 GeV, respectively. The gaugino-to-scalar mass ratio is r g = 3.5. 

Note that h\ is an (9(10 -3 ) quantity as assumed above. The resulting abundance of gravitinos 
produced at the time of reheating is given by 

n$ * * 0.016 x ( ^) 1/3 r^r 6 r m * ) ~ 1/2 7 (^-) ~ 4 , (4.i9) 

3/2 \9° 8C J ^3.5/ V300GeV/ Vl.5TeV/ V15GeV/ ' V ; 

which is required to be below the present dark matter abundance. This gives further constraint 
on the parameters in our model. 

Now, we see the existence of the parameter space where all the constraints are satisfied. 
Figure [3] shows the allowed region (shaded region) of the gravitino mass (horizontal axis) and 
the gluino mass (longitudinal axis) which satisfies the constraints from the decay temperature 
f)4.12p and the gravitino abundance (I4.17|) . (14.191) . We also show the parameter region where 
the total gravitino abundance corresponds to the present dark matter abundance. Here, the 
moduli mass is taken to be 300 GeV (left figure) or 500 GeV (right figure). The gaugino-to- 
scalar mass ratio is set to be r g = 3.5 in both cases. The constraint 2h% > hy is presented so 
that the inflaton and the waterfall fields would not decay into the messengers, which would 
change our thermal history. We can see from the left figure that the thermally produced 
gravitino abundance and the decay temperature put strong constraints on the gluino mass. For 
the region of the small gravitino mass, the abundance of non-thermally produced gravitinos 
has an important effect on the allowed gluino mass region. On the other hand, in the right 
figure, the abundance of gravitinos from the pseudomoduli decay puts a strong constraint on 
the allowed region of the gluino mass and the constraint from the decay temperature becomes 
weak. The allowed gravitino mass is always larger than O(10) GeV. It may be possible for 
the gravitino mass to be much larger than O(10) GeV, but it is not favored because the 
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5 10 15 20 25 30 20 40 60 80 

m 3 / 2 /GeV m 3/2 /GeV 

Figure 4: The allowed region (shaded region) of the gravitino mass (horizontal axis) and the 
gaugino-to-scalar mass ratio (longitudinal axis). In the left and right figures, the moduli mass 
is taken to be 300 GeV and 500 GeV respectively. The gluino mass is 1.5 TeV. 

flavor-changing neutral current would not be suppressed in such a case. Then, we conclude 
that the gravitino mass of the order of 10 GeV is favored for our model. We define the ratio 
of the gravitino abundance produced at the time of reheating and the pseudomoduli decay, 
r 3 / 2 = ^3/2* /^3/2 - ^ S ^ ne S m i no niass is decreasing, the ratio also decreases and the abundance 
of gravitinos from the pseudomoduli decay becomes dominant. 

Figure H] shows the allowed region (shaded region) of the gravitino mass (horizontal axis) 
and the gaugino-to-scalar mass ratio (longitudinal axis). In the left (right) figure, the moduli 
mass is taken to be 300 GeV (500 GeV). Here we take the gluino mass as 1.5 TeV. In the 
left figure, the upper bound for the mass ratio comes from the constraint that the moduli 
oscillation is stable, in other words, the messenger fields do not become tachyonic during 
the moduli oscillation. The constraint from the gravitino abundance from the pseudomoduli 
decay gives a weaker bound. On the other hand, in the right figure, the gravitino abundance 
from the pseudomoduli decay gives the stronger constraint on the mass ratio. We can also 
see that the abundance of gravitinos produced at the time of reheating gives a lower bound 
for the mass ratio in both the left and right figures. As the mass ratio is decreasing, the ratio 
r3/2 increases and the abundance of thermally produced gravitinos becomes dominant. 

As discussed in Ref. [20], if the Bino is the next to lightest supersymmetric particle (NLSP), 
the moduli field can decay into two Binos. In this case, these Binos decay into gravitinos later, 
which breaks BBN. However, in the case that the moduli, gluino and gravitino masses take 
the reference values and the gaugino-to-scalar mass ratio is r g = 3.5, the Bino mass is around 
230 GeV. Then, the moduli field cannot decay into two Binos and hence this problem does not 
occur. In conclusion, there is a model parameter space where this scenario can be successful. 
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5 Conclusion 



In this paper, we have proposed an inflationary universe scenario with gauge-mediated SUSY 
breaking. The higher mass scale in the model is set for the inflationary scale, and the lower 
mass scale corresponds to the SUSY breaking scale to give the correct MSSM soft masses 
by direct gauge mediation. After inflation, the metastable SUSY breaking vacuum is chosen 
naturally. We have analyzed the reheating stage of the model. We have studied the pseudo- 
moduli oscillation and its decay as was studied in Ref. [20]. We find that the pseudomoduli 
oscillation dominates the universe after reheating and dilute thermally produced gravitinos. 
Gravitinos are also produced by the pseudomoduli decay. We find a model parameter space 
consistent with the gravitino dark matter. 

We would like to comment on the two issues remained in our model. One is the cosmic 
string problem. We here take a hybrid inflation model as a concrete realization of our cos- 
mological scenario, which suffers from the cosmic string problem. We have set the coupling 
constant hy as C>(10~ 3 ) in order to suppress the gravitino abundance produced directly from 
the inflaton decay. Then, the cosmic strings may be produced at the end of inflation be- 
cause the waterfall fields break the C/ (l)i symmetry and their tension G[i is calculated as 
Gfi ~ lum} jhy. Recent observation of the CMB constrains the value of the cosmic string 
tension, Gfi < (2 — 7) x 10~ 7 [28]. In the case of the standard hybrid inflation with the loga- 
rithmic potential, the constraint on the coupling constant is expressed as hy <C 1CT 5 , which 
is conflicted with the typical value we have applied^ Therefore, the cosmic string problem is 



severe in our present model. However, we can modify our model to avoid this problem [30] by 
changing the vacuum structure [31] or replace the model from the standard hybrid inflation 
model to the shifted or smooth hybrid inflation model [32] • By applying them, the cosmic 
string would become unstable or be diluted enough not to affect the CMB or GW observation. 
We leave the categorization of these modifications for a future work. 

Another issue is baryon asymmetry in the present universe in our scenario. Since there is 
a large dilution factor A -1 ~ 10 -3 , we have to produce a large amount of baryon asymmetry 
before the pseudomoduli domination. One candidate is the Affleck-Dine mechanism [33jo It 
may be interesting to find a way to generate baryon asymmetry in the SUSY breaking sector, 
which is also left for the future study. 



12 Another constraints on the cosmic string tension from the gravitational waves (GWs) has been reported 
recently [29] . 

13 The suppression of the net baryon asymmetry due to the nonzero temperature effect |34j can be avoided 
if we consider the multiple field case [35] . 
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Appendix 

We can also consider the case where the moduli mass is smaller than the twice the Higgs 
mass, m$ < 2rrih, and the moduli field decays dominantly into a gluon pair. The interaction 
with the standard model gauge boson can be extracted from the moduli dependence of the 
running gauge coupling #3(3>). The result is given by 

1 gl hyhq, mz 
4 (47r) 2 h\ m? 

Then, the decay rate of this process is calculated to be 



As in the case discussed in the main text, the temperature after the moduli decay can be 
calculated from the above decay rate as 

T d * - J= V2 ™^ 2.3 MeV x (^L_f 2 (J*L.\ (A . 3) 

V M Pi/ ™3/2 VIOGeV/ VlTeV/VlMeV/ V ; 

The temperature is required to be above 2 MeV so that the standard BBN properly occurs. 
This requirement constrains on the parameters in our model just as in the case in the main 
text. By the way, using the partial decay width (I4.14p . the branching fraction of the decay 
into gravitinos is given by 

Then, the present energy density-to-entropy ratio of gravitinos is estimated as 
P3/2 _ 1/3 m„\ . . in _ 12 _, / ™<f \ 5 / 2 / rn- 9 
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^4.4xlO- 12 GeVxf-^-) 5/2 (-^-)" (A.5) 
\M Pl J m~ g VIOGeV/ VlTeV/ v ; 



22 



On the other hand, although gravitinos produced at the time of reheating is diluted by the 
pseudomoduli decay, their present energy density-to-entropy ratio is still large, 

Thus, from (I3.23p . we can see that the total present energy density-to-entropy ratio of grav- 
itinos, + P3/2 V s ' * s so l ar g e that they overdose the universe. As a result, this scenario 
cannot be compatible with the present universe. 
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